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SUMMARY 


A free-f light investigation was conducted to determine the aerodynamic char- 
acteristics and motions of a rocket-boosted model of a simplified hypersonic 
glider configuration with wing-tip fins and a blunt base through a Mach number 
range extending up to 4.4. The model was statically and dynamically stable at 
small angles of attack and experienced coupled motions when pitched to an angle 
of attack of about 12°. The addition of the wing-tip fins resulted in no improve- 
i ment of the vehicle motions compared with the same vehicle with no fins. 

In general, the coupled motions experienced by the model were simulated by 
the use of the equations of motion for five degrees of freedom. The analog study 
showed that the rate of change of effective-dihedral derivative with change in 
angle of attack was a significant cross- coupling derivative for the configuration 
studied and should be considered for flight prediction of the motions for config- 
urations which are highly swept and have the mass concentrated along the longitu- 
dinal axis. An analytical solution for a simplified three-degree-of-freedom set 
of equations yielded results that compared satisfactorily with the flight records 
of rolling velocity and yielded a rolling velocity parameter which affected the 
type of rolling motion experienced. The appendix presents the analytical solution. 


INTRODUCTION 


The current interest in highly swept, delta wing configurations for high- 
speed reentry and manned orbital flight has prompted a free-flight investigation 
to determine the aerodynamic characteristics and dynamic behavior of several- 
research models with these general characteristics. Reference 1 presents the 
results of one such test and shows that the model angular velocities became large 
and erratic following a disturbance. The rolling velocities were of particular 
interest, with periodic positive and negative values along with short-period 
sinusoidal-like peak-magnitude variations of 35 radians per second being exper- 
ienced. Electronic analog motion simulation showed that the primary reason for 


this behavior was the concentration of mass along the axis of symmetry, the low 
roll inertia, and the large increase in rolling moment due to sideslip with 
increasing angle of attack. 

The object of this investigation is to present the dynamic behavior of a 
delta wing configuration which is the same as that of reference 1, except for the 
addition of wing-tip fins which were added with the idea of increasing directional 
stability, of reducing coupling forces, and hence, it was thought, of decreasing 
the roll rates and accelerations. 

The model was flight tested at the NASA Wallops Station over a Mach number 
range from 4.4 to approximately 1.5 with corresponding Reynolds numbers, based 
on wing mean aerodynamic chord, of 42.5 X 10^ and 0*7 X 1C)6. Pulse rockets were 
used to disturb the model in pitch to an angle of attack of about 12°, and a 
10- channel telemeter was used to transmit continuous acceleration and pressure 
data to ground receiving stations. 

Aerodynamic coefficients based on accelerations of the center of gravity are 
presented. With the use of f ive-degree-of-f reedom equations of motion, an analysis 
of the rolling motions was made by comparing the free-flight data with analog runs 
for similar altitude and aerodynamic- coefficient conditions. In addition, an 
analytical solution of a simplified set of three-degree-of-f reedom equations of 
motion were used to predict rolling motions. 


SYMBOLS 


The basic data in this report are presented with respect to an axis system 
originating at the 54.3-percent station on the longitudinal axis of the rocket 
model. (See fig. 1.) 

a T longitudinal accelerometer reading referenced to center of gravity, 

g units 


%Hi 


high- range longitudinal accelerometer reading, positive in positive 
x-direction, g units 


a L,Lo 


low-range longitudinal accelerometer reading, positive in positive 
x-direction, g units 


%cg 


normal accelerometer reading referenced to center of gravity, g units 


a N n normal accelerometer reading at nose location, positive in negative 

z-direction, g units 

a^ ^ normal accelerometer reading at tail location, positive in negative 

’ z-direction, g units 


a y 


lateral accelerometer reading, positive in positive y-direction, g units 
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^Cg 

b 

CD 


C I<x 


' 3 a 
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' m a 


■un. 

a 


lateral accelerometer reading referenced to center of gravity, g units 


wing span, 1.77 ft 
drag coefficient. 


Drag 

gas- 


lift coefficient, Cp « Cjj cos a, 


Lift 


lift- coefficient derivative, 
rolling-moment coefficient, 
damping-in-roll derivative. 


per radian 
da 

Rolling moment 


dCi 




OcoSb 


, per radian 


rate of change of rolling-moment coefficient with change in yawing 


angular velocity factor. 


dC, 




2 V 


■, per radian 


dCn 


effective-dihedral derivative, — -, per radian 

rate of change of effective-dihedral derivative with change in angle 
^ c Zr 

of attack, per square radian 

da 

pitching-moment coefficient, Pitching^momen t 

^00 

rate of change of pitching-moment coefficient with pitching angular 


dC- 


velocity parameter, — 7— Er P er radian 

m 

dC 

static stability derivative, — per radian 

da 

rate of change of pitching-moment coefficient with rate of change of 

dCrri 

angle-of-attack parameter, — /. v -, per radian 
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Cmr: 


rate of change of pitching-moment coefficient with angle of sideslip. 


be 

W 


m 


, per radian 


C N 

C N a 

Cn 


normal-force coefficient. 


Normal force 




Sc 


N 


normal-force-coefficient derivative, — -, per radian 

oa 

. . . Yawing moment 

yawing -moment coefficient, — 

q^Sh 


Cn r 


rate of change of yawing-moment coefficient with rolling angular 

Sc r 


velocity parameter. 


"n 


bl^ 


2V 


, per radian 


Cn T 




rate of change of yawing-moment coefficient with yawing angular 

sc. 


velocity parameter. 


n 


sAk 

\2V 


, per radian 


Sc 


directional- stability derivative, per rac ^- an 


"n- 


rate of change of yawing-moment coefficient with rate of change of 


Sc. 


n 


angle-of-sideslip parameter, P er radian 

s& 

\2V 

q^Sb p 

° lx ° 

Cv axial-force coefficient, positive in positive x-direction. 


Axial force 


lateral-force coefficient. 


Lateral force 




c y t 


rate of change of side-force coefficient with angular velocity factor 

SCv 

in yaw, — 7 \ , per radian 

SCv 

lateral-force-coefficient derivative, , per radian 

oP 


k 


c 

F,K 

g 

h 


ix 

Iy 


Iz 

k 

l 

M 

m 

P 

Pb 

Po 

Poo 

<1 

<3oo 

R 

r 

S 

t 

V 

W 

X,Y,Z 


wing mean aerodynamic chord, 3.83 ft 

incomplete and complete elliptic functions of first kind 
acceleration due to gravity, 32.2 ft/sec 2 
model altitude, ft 

mass moment of inertia in roll, ( slugs )(sq ft) 
mass moment of inertia in pitch, ( slugs )(sq ft) 
mass moment of inertia in yaw, ( slugs )(sq ft) 

critical roll- rate parameter (eqs. (Al 4 ) and (A 23 )) 
reference length, ft 
Mach number 

model reference mass, w/g, slugs 
rolling angular velocity, radians/sec 
base pressure, lb/sq in. 

total pressure, measured at nose stagnation point, lb/sq in. 


free-stream pressure, lb/sq in. 

pitching angular velocity, radians/ 

pV 2 

free-stream dynamic pressure, 


sec 

lb/sq ft 


Reynolds number, based on wing mean aerodynamic chord ( 3*83 ft), 


pVZ 

P 


yawing angular velocity, radians/sec 

wing reference area, 5*^7 sq ft 

time, sec 

velocity, ft /sec 

model weight, lb 

fixed body- axis system 

distances measured along the X-, Y-, and Z-axes 
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a angle of attack, radians 

0 angle of sideslip, radians 

7 flight-path angle referenced to local horizontal, deg 

[i free- stream coefficient of viscosity, lb-sec/sq ft 

p free-stream density, slugs/cu ft 



0 angle of roll, radians 

\| f = sin" 1 (k sin 0) 

Subscripts: 

min minimum value 

max maximum value 

o initial condition 

A dot ( * ) over a symbol denotes a first derivative with respect to time, and 
double dots (") over a symbol denote a second derivative with respect to time. 

A bar (") indicates an angle normalized with respect to o Q . 

MODEL AND INSTRUMENTATION 


Model 

A dimensional drawing of the rocket-boosted model is shown in figure 2 . The 
reference areas and lengths and physical characteristics of the model are given 
in table I. Photographs of the model are presented in figure 3 j and a photograph 
of the model and booster in launch position is shown in figure k. 

The model had planar symmetry in the horizontal and vertical planes, sharp 
leading edges, and a blunt trailing edge. The planform of the basic body was a 
78.87° sweptback clipped delta surface with streamwise tips, aspect ratio of 0.53* 
and a taper ratio of 0 . 191 . The wing-tip fins were clipped delta surfaces swept 
back 70 °. The vertical tail had a sweepback angle of 76°. The planar surfaces 
of the model were not properly alined in some instances because of tolerances of 
construction giving rise to slight aerodynamic asymmetries. A detailed descrip- 
tion is presented in table I and figure 2 of the physical and mass characteristics 
of the complete model and of each component. For all practical purposes, the mass 
balance of the model was such that the principal axes were coincident with the 
body axes of symmetry. 
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Instrumentation 


The rocket-boosted model was equipped with a 10-channel telemeter which con- 
tinuously transmitted information concerning longitudinal, transverse, and normal 
acceleration; angular acceleration in pitch and yaw; angular velocity in roll; 
and total and base pressures. There were two longitudinal accelerometers, one 
with a g-unit range from 1 to -12 and one with a g-unit range from 1 to -2. There 
were also two normal accelerometers, one mounted in the nose section of the model, 
and the other mounted in the rear of the model. The base pressure data repre- 
sented an average over the semispan of the blunt trailing edge of the model as 
obtained by a manifolded tube, as shown in figures 2 and 3(b)* The locations of 
the instruments with respect to the center of gravity are given in table II. 

The data accuracy of the 10 instruments installed in the rocket-propelled 
model was approximately ±5 percent of the calibrated scale of each instrument. 

The error in the aerodynamic coefficients was dependent on the free-stream dynamic 
pressure, which had a possible variation of 25 percent, on velocity, which had 
possible variations of ±100 feet per second, and on density, which had possible 
| variations of ±5 percent. 


TESTS 


Figure 4 is a photograph of the rocket model and launch vehicle elevated on 
i the launcher prior to firing at the NASA Wallops Station. The three- stage launch 
vehicle consisted of an Honest John rocket with four standard 12- square-foot 
trapezoidal fins for a ground-fired first stage; a Nike M5 rocket motor with four 
standard 5-square-foot trapezoidal fins for a delayed, ground-fired second stage; 
and a Nike with four standard 5- square-foot 4° wedge fins for a ground initiated, 
timer-fired third stage. 

A small pancake rocket motor was installed in the base of the model behind 
the blowout diaphragm (see fig. 3 (b)) and was programed to separate the model 
from the burned-out third- stage Nike. Preflight calculations of the relative 
| decelerations of the burned-out third-stage Nike and the model, based on drag- 
weight ratios, indicated the need for an additional separation force in order to 
prevent collision after separation. In order that this requirement be satisfied, 
two 3.25-inch Mk 7 rocket motors were strapped in reverse thrust position to the 
forward end of the third-stage Nike, as shown in figure 4. The nozzles of the 
two rocket motors were canted out 20° with respect to the center line of the third 
stage and were plugged and sealed to protect the grain and igniters from aero- 
dynamic heating during the boosted portion of the flight. 

Two pulse rockets with 27 pound-seconds of impulse (burning time 0.05 second) 
were located in the nose section of the model (see fig. 2) and so oriented that 
the direction of thrust was in the positive and negative z-directions. Approxi- 
I mate calculations prior to the free-flight test for which two- degree- of- freedom 
pitch-plane dynamics and aerodynamic coefficients from reference 2 were used 
indicated that the maximum disturbance in the pitch plane alone would be about 
12 °. 
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Continuous flight-path data were obtained by tracking the model with three 
radar sets: an AN/FPS-16, an SCR-584, and an SCR-584 Model II. Atmospheric con- 

ditions were obtained from a Rawinsonde balloon which was released prior to launch. 
Mach number and dynamic pressure were determined from the total pressure measure- 
ments of the model and the ambient pressures and temperatures at corresponding 
altitudes along the flight path. The error of the most accurate tracking radar 
set used ( AN/FPS-16) was 0.1 mil for the given angular measurement and 50 feet 
for the given linear measurement. 

Data are presented in this report ranging from a Mach number of about 4.4 at 
an altitude of about 34,000 feet to a Mach number of 1.1 at an altitude of about 
108,000 feet. The variation of altitude with Mach number is presented in fig- 
ure 5 . Variations of dynamic pressure and Reynolds number with Mach number are 
presented in figures 6 and 7 > respectively. 

Plots of 10 channels of telemeter data, dynamic pressure, and Mach number 
for three time intervals of the flight test of the model are presented in fig- 
ure 8. The points were machine plotted from magnetic tape records and show 
occasional scatter points due to noise. Figure 9 presents an envelope of the 
variation of maximum roll rates plotted against Mach number for the flight. 


ANALYSIS 


Basic Data Reduction 

The measured free-flight linear-acceleration data were corrected to the cen- 
ter of gravity because linear accelerometers which are not mounted on the center 
of gravity measure not only the translatory accelerations but also the accelera- 
tions due to angular velocities and angular accelerations. The following 
relationships were used to correct the accelerometer readings to the center of 
gravity: 


%cg - %Lo + g 


%Lo 


(l 2 + r 2 ) + ya L)Lo (r - Pi) - z a T ,. Il0 ( 4 + P r > 


'L,Lo 


( 1 ) 


a Y,cg = ay + 


yay( r2 + p2 ) + z a Y ( p " < l r ) " x a v ^ + pt l) 


V ay v 


( 2 ) 


and 


a N, eg " % n + g 


" Za N,n( 


P 2 + <1 2 ) - 


Xa N,n (4 " pr) + y a N , n (i> + qr) 


(3) 


where the distances x, y, and z from the center of gravity to the individual 
accelerometer are given in table II. The variations of r, and p were meas- 

ured directly from the angular accelerometers and the roll gyro, respectively, 
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and required no corrections. The variations of q and r were obtained by- 
integration of 4 and r over several cycles and by assuming that they varied 
symmetrically about zero values. The angular acceleration in roll p was 
obtained by differentiating p with time. 

The values of a L Lo and ajj n were used for determining a^ cg and a^ cg , 
respectively, throughout most of the flight. The values of a L^Hi were used only 
at the higher Mach numbers where the deceleration was greater than the range of 
a L Lo ; and a N t served only as a check for a f^ cg * 

The total force coefficients were determined from the instantaneous values 
of the translatory accelerations (referenced to the center of gravity) as follows; 

(*0 


(5) 


(6) 


w 


C N - %cg 


_W_ 


" C x - — a L, cg 


and 


W 


C Y " 4^s a Y, eg 


The total moment coefficients for the rocket model were computed from the instan- 
taneous values of angular acceleration and velocity by using the following 
expressions: 


Cm 




(7) 


and 


C n 


-^-r 


+ 


fcy - x x\ 
V ) 


PT 


( 8 ) 


c z 


i2Lp 

q~sb p 


+ 


Iz - ^ 

< J«, Sb / 


(9) 
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Analog Simulation 


An analog simulation analysis similar to that of references 1 and 3 was used 
to study rolling velocities of the model near Mach numbers 4 and 2. There were 
three reasons for using the analog computer for motion simulation: (l) the vehicle 

aerodynamics were known to be nonlinear with angle of attack or yaw; (2) the 
lateral and longitudinal aerodynamics were coupled by the rolling velocity; and, 

(3) the vehicle motions were not symmetrical. The equation associated with the 
drag force was omitted since the velocity of the model was nearly constant through 
the time intervals studied. The cross-product inertia terms have been omitted 
from the equations because the principal axes coincided with the body axes. The 
equations of motion were written for the body-axis system and are as follows: 

Normal force: 


i = 1 - PP - ^ C^a 


(10) 


Side force: 


* = " r + ap + ST ‘V 


(11) 


Rolling moment: 


P = 


~ I Z\ , q oo Sb 


qr + 1 (c lp (> ♦ ♦ |jfp(c !p P + C Ir r) (12) 


Ix 


Pitching moment: 


4 = f z i Y x ) pr + Cm a a + Mr( Cm( i q + Cm&d 


(13) 


Yawing moment: 


r = 


d x - i Y 




P<1 + c n R P + ^T-( c n r r + C n -£ + c np p) 


IZ I Z n P H 2VI Z 


(14) 


The body was assumed rigid and the gravity terms were omitted from the equations. 
Certain aerodynamic terms, such as C m o and Cy r ^ were omitted from equations (10) 
to (l4), inasmuch as their effect on the motions was considered to be negligible. 

In the analog study it was assumed that p = -r in equation (l4), thus allowing 
some simplification in the last term. 

The equations of motion for five degrees of freedom were used at M = 4 and 
M = 2. The static stability derivatives from reference 5 were used in the analog 
study at M = 2 after they were transferred from the reference point of refer- 
ence 5 (66.6 percent of the theoretical body length behind the theoretical apex) 
to the center of gravity of the vehicle. The damping- in-pitch coefficients Cm^ 
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and were initially calculated from linearized theory (refs. 6 and r j), but 

their values were lowered when they were found to be excessive for an accurate 
simulation of the motions. 


The damping- in- yaw coefficients Cn r and Cnjj were also calculated in a 

manner similar to that for pitch, with the assumption that all damping was due to 
the wing- tip fins and vertical tail. The value for C z was estimated from ref- 

erence 8, and values for C 7 and C n were estimated from lifting-line theory. 

Table III presents values for all the derivatives used. It should be noted that 
two slopes were used to approximate Cyp, C n p, and because they were non- 

linear with p and/or a. (These slopes were obtained from wind-tunnel tests of 
ref. 5.) The initial conditions for the equations of motion were obtained at a 
given time in the free-f light record. Initial values of a and P were found 
by solving the equations of motion. The initial conditions for the analog inves- 
tigations made for M = 4 and M = 2 are given in table IV. 


Because of the large predominating effect of C Zf) on the motions, a highly 

Pa 

simplified motion study was made with the use of the following three-degree-of- 
freedom equations of motion: 


d = -Pp 

(15) 

UD* 

II 

£ 

( 16 ) 

q^sb 

-f — C Zr pa 
I X Pa 

(17) 


An analytical solution for this simplified set of equations is included in the 
appendix . 


RESULTS AND DISCUSSION 


Basic Data 


One of the purposes of the flight test presented in this paper was to deter- 
mine if the roll rates and accelerations of the basic vehicle could be decreased 
by the addition of two wing-tip fins. Reference 5 presents longitudinal and 
lateral stability data for the model of reference 1 (the basic model) at Mach num- 
bers of 1.4l and 2.01 and for the model of this report at a Mach number of 2.01. 
The Reynolds numbers of these wind-tunnel tests are given in figure 7- As shown 
in reference 5* the vehicle with wing-tip fins was statically stable at a Mach 
number of 2.01. As a result of adding wing-tip fins to the basic configuration. 


wind-tunnel data (ref. 4) at M = 2.01 
from 0.0005 to 0.0022 per degree. 


Pa 


and a = 4° show that C n p was increased 
was decreased from -0.0012 to -0.0014, 
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CSfp was decreased from -0.0020 to -0.0045, and C m was decreased from 0.0025 

to - 0 . 0010 , all of which were referenced to the 66.6 percent longitudinal station 
behind the theoretical apex and/or the wing mean aerodynamic chord. Even though 
the directional and longitudinal stability was increased, a comparison of the 
roll-rate histories for the present test and those of reference 1 shows similar 
types of oscillations for both models. The reason for this result is that the 
addition of the wing-tip fins did not result in a decrease in ^Zp a # 

Figure 8 (a) shows the time interval from 23.5 to 24.5 seconds, which begins 
just after model separation (23*48 seconds) and includes the period during which 
the first pulse-rocket ignition occurred ( 23.52 seconds). This figure indicates 
that the model experienced oscillations of large amplitude, due to the pulse 
rocket, which were coupled in pitch and yaw. The roll rate oscillated about zero 
through large amplitudes during and immediately after pulse-rocket ignition and 
then steadily increased to a large positive amplitude of approximately 35 radians 
per second. The second pulse rocket was ignited at 27 . 01 seconds. Figure 8 (b) 
presents data from 28.5 to 29-5 seconds and shows the sinusoidal oscillation of 
the model in roll, pitch, and yaw. Figure 8 (c) presents data from 38.6 to 
39.6 seconds for M = 2 after ignition of the second pulse rocket. In general, 
the motions of the model were coupled, and the roll rate over some intervals was 
unsteady and varied randomly from positive roll oscillations through roll 
reversals . 


Vehicle oscillation about the stability boundary was considered as a pos- 
sible cause of the coupled motions. By using the methods of references 4 and 9, 
it is possible to superimpose a region of roll divergence for steady roll rates 
on the curve of figure 9 . The figure indicates that from M=4.2 to M = 3-7 
the roll rates are not in the regions of divergence; however, at lower Mach num- 
bers the roll rates cross back and forth over the boundaries of divergence. It 
should be noted that the moments are coupled prior to encountering the regions of 
divergence, which indicates that the cause of the coupling is purely aerodynamic. 
It should also be noted that the boundaries were derived for constant roll rates 
and linear equations of motion, which do not account for the second-order effects 

such as that of C 7 Q , and, consequently, the boundaries may not be applicable for 

Pa 

this flight condition. 


The variation of the total force and moment coefficients with time is pre- 
sented in figures 10 to 13 for four time increments. The time increments were 
picked to show variation in aerodynamics (l) after model separation and during 
first pulse-rocket ignition at M 4.1; (2) after first pulse-rocket disturbance 
at M ^ 3.9; ( 3 ) after second pulse-rocket disturbance at M = J.l; and (4) after 
second pulse-rocket disturbance at M « 2.1. Parts (a) and (b) of these figures 
are the longitudinal coefficients, total normal force and total pitching moment, 
respectively. Parts (e) and (f) are the directional coefficients, total side force 
and total yawing moment, respectively; parts (c) and (g) are total axial force and 
total rolling -moment coefficients, respectively. Parts (d) and (h) of figures 10 
to 13 are variations of total pitching moment plotted against total normal-force 
coefficient and of total yawing moment plotted against total side-force coeffi- 
cient, respectively. A line representing the approximate slope of the data has 
been drawn through some of the curves indicating static stability as shown by 
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reference 5 but these curves do not pass through zero (at Cjj = 0 and Cy - 0), 

probably because of the asymmetries due to model construction, as well as the 
effects of the cross- coupling moments and forces. The model also appears to 
exhibit a trim change with Mach number. 

In figure 10(c) the reduction in axial force is evident and is due to after- 
burning in the separation motor from about 23-62 to 23*90 seconds. The variation 
in base pressure shown in figure 8(a) is also due to the effects of afterburning. 


Analog Analysis 

Five-degree- of -freedom analog analysis .- The first attempt at analyzing the 
flight data herein was to compare the longitudinal and lateral aerodynamic coeffi- 
cients and derivatives obtained by simple analysis with those presented in refer- 
ences 1 and 5. However, the coupled forces and moments did not allow sufficient 
steady- state motion for comparison. 

Numerous analog runs were made using the f ive-degree-of-freedom equations at 
M = 4 and M = 2. In an effort to make adjustments to the coefficients for 
proper motion simulation, it was found that c Zp a was predominant in its effect 

on the rolling motion. By the use of the aerodynamic coefficients of table III 
and the initial conditions of table IV, figures l4 and 15 were selected as repre- 
sentative of the flight rolling velocity at M = 4 and M = 2, respectively. 
Figures l6 and 17 are reproductions from the flight record of the rolling velocity 
and pitch and yaw accelerations for M = 4 and M = 2, respectively. 

Comparison of figures 14 and 16 indicates that a qualitative type of agree- 
ment is obtained between the analog record and the flight record. Comparison of 
figures 15 and 17 indicates the same type of agreement, but in this case there 
appears to be closer quantitative agreement. 

Three-degree- of -freedom analog analysis .- The manner in which Cjp^ seemed 

to dominate the five -degree- of -freedom analog studies indicated that a greatly 
simplified approach to the simulation of the rolling velocities was possible. 
Equations (15), (l6), and (17) were used for the simplified three-degree-of- 
freedom analysis . Analog studies were made of the rolling velocities at M = 4 
(fig. 18) and at M = 2 (fig. 19). Comparison of figures l8 and 14 shows the 
effect of the simplified analysis. In figure 18 the frequency and amplitude of 
the rolling- velocity oscillation is more obvious than in figure l4 and agrees 
better with the periodic rolling velocity shown in the flight data at M = 4 
(fig. 16). Figure 19 is an analog record of the rolling velocity at M = 2 and 
shows, in the initial moments of the run, a roll-reversal motion similar to those 
of the flight record followed by a negative periodic rolling velocity. Figure 19 
was selected as being typical of the analog runs made at M = 2 with the simpli- 
fied equations of motion. 

Consecutive analog runs showed that the simulated motions would not repeat 
exactly for identical initial conditions. The angular velocities and their fre- 
quencies were similar to those from the test. Small variations in the stability 
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derivatives and use of other sets of initial conditions did not provide repeat- 
ability. The nonrepeatability appeared to be due to the sensitivity of the equa- 
tions to small computer errors of the analog for the inputs used. These small 
computer errors may be compared to spurious disturbances during the flight and, 
in this case, are partly responsible for the realistic simulation of the free- 
flight motions. The analytical solution to the equations of motion given in the 
appendix shows that the type of rolling motion experienced depends on a param- 
eter k which in turn depends on the initial conditions, dynamic pressure, and 

the Ci n derivative. For values of k < 1, the solution gives rise to a peri- 
</ pa 

odic positive or negative rolling velocity, whereas, for values of k > 1, oscil- 
latory rolling velocities (roll reversals) are predicted. A value of k = 1 
yields a logarithmic decay to zero of rolling velocity, but it is a very sensitive 
condition and not likely to show itself in flight. 

Several calculations have been made using the analytical expressions of the 
appendix to make comparisons with flight rolling- velocity time histories. The 
rolling- velocity time histories for k < 1, k = 1, and k > 1 were computed for 
the comparisons. The value of k for the case in which k < 1 was determined 
from equation (A23) of the appendix for which the flight values of p ma x and - 
Pmin were used; hence, in this instance, the three-degree-of-freedom equations 

of motion may be judged only by the agreement of the general shapes of the flight 
and computed curves and the period of the motions. Figure 20 (k = 0.80) shows 
that the simplified three-degree-of-freedom equations predicted both the period 
and general shape of the rolling- velocity variations very well and also indicates, 
at least for the portion of the flight shown in figure 20, that the q and r 
velocities did not appreciably affect the rolling motion. 

Calculations of k for k = 1.004 were made with the expression for the 
period of oscillation given as equation (A2l) in the appendix. Figure 21 shows 
that the amplitudes are fairly well predicted by the three-degree-of-freedom 
equations of motion; however, the shapes of the two curves do not agree as well 
as for the case when k = 0.80. 

A portion of the flight record which has a value of k of approximately 1.0 
has been selected for comparison with the theoretical prediction in figure 22. 

For the actual theoretical case, the half- cycle logarithmic decay of roll rate to 
zero occurs at infinity. For practical purposes, two half-cycles which decay 
logarithmically to zero and which have a half-period equal to flight-record half- 
period obtained for k = 1.0 from the three-degree-of-freedom equations of motion 
are compared with the flight record. Figure 22 shows that the comparison can be 
considered comparable in only a qualitative manner. The differences between the 
theoretical curve and experimental values may be due largely to the fact that the 
simplified theory neglects the effects of q and r. Better agreement is indi- 
cated in figure 15 where the rolling-velocity time histories have been determined 
from the equations of motion for five degrees of freedom and for values of k of 
approximately 1.0. 

It should be noted that for similar highly swept low-aspect-ratio glider 
configurations with and without planar symmetry, as described in the wind-tunnel 
tests of reference 5, is approximately the same magnitude. According to 

Pa 
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the general analog study made herein, the aforementioned configurations would he 
subject to similar coupled motions as obtained from the present flight test at 
supersonic speeds unless provision is made for roll or yaw control. 


CONCLUDING REMARKS 


Results from a free- flight investigation conducted to determine the aerody- 
namic characteristics and motions of a simplified hypersonic glider configuration, 
which had planar symmetry and the mass concentrated along the longitudinal axis, 
indicated that the flight motions were coupled and nonlinear and that the model 
experienced unsteady rolling Velocities up to 35 radians per second. The direc- 
tional and longitudinal stability was increased by the addition of wing-tip fins, 
as shown by wind-tunnel tests. The coupled rolling motions of the model were 
similar to those of the same model without the wing- tip fins. The model was stat- 
ically and dynamically stable in pitch and yaw at small angles of attack and side- 
slip up to a Mach number of 4. 

Electronic analog studies using equations of motion for five degrees of free- 
dom of the rolling velocities of the model indicated that the rate of change of 

the effective- dihedral derivative with change in angle of attack C, was the 

Pa 

predominant term in the simulation of the amplitude and period of the rolling 
motions. An analytical solution of the simplified three- degree-of-freedom equa- 
tions of motion containing the derivative yielded a critical roll- rate 

Pa 

parameter which correlated the type of rolling motions experienced by the 
configuration. 

The analog study showed that the term C 7 was a significant cross- coupling 

Pa 

derivative for the configuration studied and should be considered for predicting 
the flight motions of configurations which are highly swept and have their mass 
concentrated along the longitudinal axis. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va. , July 30 , 19 ^ 3 • 
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APPENDIX 


AN ANALYTICAL SOLUTION OF THE SIMPLIFIED 
TKREE-DEGREE-OF-FREEDOM EQUATIONS 
OF MOTION 

By Percy J. Bobbitt 


The simplified three-degree-of-freedom equations 


a = -Up = -P f> 


(Al) 

£ = ap = a0 


(A2) 

and 

p - i x 


(A3) 

(eqs. ( 15 ), (l6), and ( 17 ), respectively, of the text) which yield motions similar 
in several respects to the flight records and five-degree-of-freedom equations are 
easily integrated and have their parallel in dynamics in the pendulum problem. Of 
primary interest here is this isolation of the critical parameter which determines 
the characteristics of the motion. 

It is a slight convenience in effecting a solution of the equations of motion 

to normalize a and P hy a 0 = y a 0 2 + P 0 2 , so that 

equations (Al) to (A 3 ) 


become 



Pi- 

ll 

S' 


(AJO 

P = ap 


(A5) 

and 



P = -0 o Pa 


(A 6) 

where 



P i» Sb p „ 2 

0 ” lx %a ° 


(A7) 

B- , ^ 

^a D 2 + p 0 2 




1 6 


and 


- a 

l/ a o 2 + Po 2 

If a is eliminated from equations (A5) and (A6), the following differential 
equation is obtained: 

pp = -CoPP (A8) 


Integrating equation (a 8) and subsequently solving for p yields 


p = 4o 2 *co?o -fcrsx 5 ) 52 


A similar equation in terms of a may be derived by eliminating p from equa- 
tions (A4) and (A6) and integrating. 


If the equation resulting from the division of equation (Ah) 
is integrated, the relation between a and p is given by 

by equation (A 5 ) 

p -p _ 2 - 2 

a + P = a o + Po = 1 

(A10) 

and it becomes clear that d and 0 may be set equal to sin 0 
respectively, or cos 0 and sin 0, respectively. Choosing 

and -cos 0, 

pi 

n 

o 

o 

CO 

(All) 

and 


3 = sin 0 

(A12) 


the equation for p (eq. (A9)) may now be written 


or 


where 



+ jl . r* t - r 

(l - kWjJ j ° 

V. ^ - 

|/Po 2 + C <A 2 


(A13) 


(A13a) 


(A14) 


IT 


The integration of equation (A13) takes three different forms, depending on whether 
k is less than, greater than, or equal to 1.0. The expressions for t are as 
follows: 


For k < 1.0, 


For k > 1.0, 


where 


and 


And, for k = 1, 


t = ± -p=|F(0,k) - F(0 o ,kj] 
yc o 


t = ± l/k) - F(i|r 0 ,l/kJ] 

yCo 


\|r = sin"-*'(k sin 0 ) 


k sin 0^1 


(A15) 


(A16) 


t = ± 



l°g e 


(l + sin 0 
\1 - sin 0 


- l°g e 


1 + sin 0o \ 

1 " sin fo) 


(A17) 


where the plus sign denotes positive values of 0 , and the minus sign denotes 
negative values of 0. (Note that F(-0,K) = F(0,k).) When 0 in equation (AI 5 ) 
becomes greater than 90 °, the extension-of -range formulas for the elliptic func- 
tion F must be used; that is, for tc/2 ^ 0 ^ Jt, 


for it ^ 0 ^ 


F(0,k) = 2F(jt/2,k) - F(jt - 0,k) 


(Al8a) 


F(0,k) = 2F(it/2,k) + F(0 - jt,k) 

and for = 0 S jt, 

F(0,k) = 4F(jt/2,k) - F(2jt - 0,k) 


(Al8b) 


(Al8c) 


For equation (Al6), f becomes equal to but does not exceed 90° ^ hence, 0 will 

not exceed sin - -*- s ^ — and 0 will reverse its sign when \|t reaches 90°; this 

condition is consistent with the fact that 0 is an even function of \|r about 
\|r = 90°. The equations for time in this case must take a different form. 
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For 


t > ± -L(F(*/2,l/k) - F(t 0 ,l/k[| 
yC o 

F(\|r,k) is replaced by 2F(jt/2,l/k) - F(t,l/k)> for 


(A19a) 


t > ± -L-|2F( rt /2,l/k) - F(* 0 ,l/k)] 
1/Co L 

F(\|r,k) is replaced by 2F(jt/2,l/k) - F(t,l/k); for 


(A19b) 


t > ± -^[^(jt^l/k) - F (\|r 0 , l/k JJ 
V^o 


(A19c) 


F(i|r,k) is replaced by 4F(jt/2,l/k) + F(\|r,l/k). 


From equations (A15) and (Al6), the revolutions per second in the case of 
k < 1 and the oscillations per second in the case of k > 1 are seen to be 


and 


Revolutions 

Second 


fa 

kkK(k) 


(k < 1) 


Oscillations 

Second 


K 

kK(l/k) 


(k > 1) 


For k = 1.0, 


p decays logarithmically to zero as 


0 goes to 90 °* 


( A20) 


(A21) 


In comparing flight roll-rate data with those computed from equations (Alj) 
to (A17), it is necessary to determine an appropriate value of k. It is possible 
to compute k directly from equation (Aik) if a number of points along a cycle 
are averaged, since 


k = 



(A22) 


However, better results will usually be obtained if k is computed from the 
formula for oscillations per second for k > 1 (eq. (A21)) or, in the case of 
k < 1, by measuring p ^y and p m -j n and noting that 


^SiS. = Ji _ k 2 or 
^max 



(A23) 
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TABLE I.- PHYSICAL CHARACTERISTICS OF ROCKET MODEL 


Wing: 

Planform area, sq ft 5*87 

Base area, sq ft 0.64 

Span, ft 1-77 

Aspect ratio 0.55 

Root chord (model length to theoretical tip), ft 5-56 

Mean aerodynamic chord, ft 5-85 

Distance from mean aerodynamic chord to root chord, ft 0.5^ 

Sweepback angle of leading edge, deg 78-87 

Dihedral, deg 0.00 

Maximum thickness ratio 0.077 

Taper ratio 0.191 

Vertical tail: 

Lateral area (two fins exposed), sq ft O .36 

Base area (two fins exposed), sq ft 0.05 

Span (total), ft O .85 

Aspect ratio 2.0 

Sweepback angle of leading edge, deg 76 . 00 

Wedge half- angle, deg 5-82 

Taper ratio 0.33 

Fins: 

Lateral area (one fin exposed), sq ft 0.40 

Base area (two fins exposed), sq ft 0.10 

Span, ft 0.62 

Aspect ratio O .96 

Planform area (two fins exposed), sq ft 0.15 

Sweepback angle of leading edge, deg 70.00 

Taper ratio 0.20 

Wedge angle, deg 7*^ 

Mass characteristics: 

Weight, lb 112.73 

Center- of- gravity position, rearward from theoretical tip, 

percent of total body length 54.31 

Mass moment of inertia in roll, (slug)(sq ft) O .38 

Mass moment of inertia in pitch, ( slugs )(sq ft) 7-86 

Mass moment of inertia in yaw, ( slugs) (sq ft) 7*99 

Product of inertia, (slug)(sq ft) 0.00 

Inclination of principal axes to model center line, deg 0.00 
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TABLE II.- MODEL INSTRUMENT LOCATIONS 


Instrument 

Quantities 

measured 

Range 

Distance from center 
of gravity 3. 




x, ft 

y, ft 

Z y ft 

Longitudinal accelerometer 

a L,Hi 

1 to 

-12 

-0.375 

0.042 

0 

Longitudinal accelerometer 

a L,Lo 

1 to 

-2 

-.375 

-.042 

0 

Normal accelerometer 

a N,n 

40 to 

-bo 

-.175 

0 

-.067 

Normal accelerometer 

a N,t 

30 to 

-30 

2.020 

0 

-.053 

Transverse accelerometer 

a y 

30 to 

-30 

-.533 

-.033 

0 

Angular accelerometer in pitch 

4 

300 to - 

300 

.050 

-.067 

0 

Angular accelerometer in yaw 

r 

300 to - 

300 

-.175 

0 

.067 

Angular velocity in roll 

P 

35 to 

-35 

.300 

0 

0 





Location 

Total pressure 

p o 

0 to 

200 

Nose tip 


Base pressure 

Pb 

0 to 

15 

Manifolded tube on 
base center line 

(fig. 3(1)) 


a Positive directions are as shown in figure 1. 
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TABLE III.- AERODYNAMIC COEFFICIENTS FOR FIVE- DEGREE- OF- FREEDOM 
EQUATIONS OF MOTION AT M = 2 AND M = k 


Coefficient 


Mach number 

M = 2 

S 

ii 

-j=- 

1.55 

0.80 

-0.20 

-0.05 

-0.01 

-0.007 

-0.01 

-0.003 

0 

0 

-0.08 

-0.005 

0.03 

0.02 

-0.05 

-0.05 

-0.05 

-0.01 


-0.15 

1 1 

0 P 

o 

-0.20 

0.28 


0.32 

-1.64 

1 1 

0 H 

H 

o o 



c i* 

Cm a 

C mq 

C >% 

C, 

Z 3 

c z 

ip 

C 7 

L r 

c np 

c n r - c np 

Cl e 

For 0 < |3 < 0.15 

For 0.15 < 3 • • 

C np 

For 0 < 3 < 0.075 

For 0.075 < 3 • • 

C Z B 

3a 

For 0 < a < 0.1 . 
For 0.1 < a . . . 



TABLE IV. 


INITIAL CONDITIONS FOR FIVE-DEGREE- OF- FREEDOM 


EQUATIONS OF MOTION AT M=2 AND M = A 


Initial conditions 

Mach number 

M = 2 

M = 4 


2.148 

5-924 

V, ft/sec 

2,096 

5,756 

Ib/sq ft 

243.7 

4,044.8 

a, radian 

- 0.196 

-0.059 

3, radian 

0.160 

-0.0022 

p, radians/sec 

13.5 

52.68 

q_, radians/sec 

-1.135 

-2.00 

r } radian/sec 

-0.10 

0 

p, radians/sec^ 

240.0 

-100.0 

4, radians/sec^ 

57-0 

222.2 

r, radians/sec^ 

21.0 

47.20 

W/qS 

0 . 08 l 

0.0049 

P*,, lL/sq in 

0.524 

2.606 

h, ft 

74,844 

42,207 




Figure 1.- Body-axis system. Axis origin is at 5 k . ^-percent longitudinal station. 

cate positive directions.) 


(Arrows indi- 
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Figure 2.- Details and dimensions of models. (All dimensions are in inches unless otherwise specified.) 



(a) Plan view of model. 
Figure 3*- Free-f light model. 


L-59-3^97 
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Figure 3.- Continued. 



Figure 3*- Concluded 
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31 


, Ib/sq ft 
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Figure 6.- Variation of free- stream dynamic pressure with Mach number. 



04 


radians/sec 2 



(a) Time interval after model separation (23A8 sec) and before, during, and after first pulse- 

rocket disturbance. 

Figure 8 .- Variations of basic data and free-stream conditions with time for three time intervals 

of flight. 
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(b) Time interval after second pulse- rocket disturbance. M = 3*0. 
Figure 8.- Continued. 
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(c) Time interval after second pulse- rocket disturbance. M = 2.0. 
Figure 8.- Concluded. 
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, radians /sec 



Figure 9.- Variation of amplitude and range of rolling velocity for free-f light model. 
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(b) Total pitching-moment coefficient. 



(c) Total axial-force coefficient. Pulse- rocket ignition occurred at 23-52 seconds. 



s 


(d) Variation of C m with Cn- 

Figure 10.- Time histories of total force and moment coefficients along with cross plots of pitch 

and yaw coefficients. M =» 4.1. 
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M > 4.1< 
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(e) Total side-force coefficient. 
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(f) Total yawing-moment coefficient. 
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(g) Total rolling-moment coefficient 



(h) Variation of Cn with Cy. 
Figure 1C.- Concluded. 
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(c) Total axial-force coefficient. 



(d) Variation of C m with 

Figure 11.- Time histories of total force and moment coefficients along with cross plots of pitch 

and yaw coefficient. M » 3 . 9 . 


bo 


- 3-92 


It » 3.84 



(h) Variation of C n with Cy. 
Figure 11.- Concluded. 


kl 


m * 3.I6 


H « 3.10 



(c) Total axial-force coefficient. 



(d) Variation of Cm with C^. 

Figure 12.- Time histories of total force and moment coefficients along with cross plots of pitch 

and yaw coefficients. M « 3 
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3.16 


* 3.W- 



(f) Total yawing-moment coefficient. 



(g) Total rolling-moment coefficient. 



(h) Variation of C n with Cy. 
Figure 12.- Concluded. 
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(a) Total normal-force coefficient. 
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(b) Total pitching-moment coefficient. 



(c) Total axial- force coefficient. 



(d) Variation of C m with C^. 

Figure 13*- Time histories of total force and moment coefficients along with cross plots of pitch 

and yaw coefficients. M « 2.1. 
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(e) Total side-force coefficient. 




(g) Total rolling-moment coefficient. 



(h) Variation of C n with Cy. 

Figure 13 . - Concluded. 
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Figure 14.- Five-degree-of-f reedom analog analysis record at M = 4. 
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Figure 15.- Five-degree-of -freedom analog analysis record at M = 2. 
4 =“ 
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Figure 1 6 .- Flight records at M = 4 of model roll rate and pitch and yaw accelerations. 
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Figure 17 .- Flight records at M = 2 of model roll rate and pitch and yaw accelerations. 




c 

ro 

•H 

T> 



Figure 18.- Simplified analog analysis record at M = 4 
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Figure 19 •- Simplified analog analysis record at 


M = 2. 
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Figure 20.- Comparison of flight data and theory for variation of rolling velocity with time, k = 0.80. 
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Figure 21.- Comparison of flight data and theory for variation of rolling velocity with time, k 


= 1.004. 
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Figure 22.- Comparison of flight data and theory for variation of rolling velocity with time, k = 1.0. 


